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Abstract 

In the context of the minimal supersymmetric extension of the Standard Model 
with the right-handed Majorana neutrinos, we study lepton flavor violating pro- 
cesses including full renormalization group running effects. We systematically com- 
pare our results with the commonly used leading logarithmic approximation, re- 
sorting to a "best fit" approach to fix all the high energy Yukawa matrices. We 
find significant deviations in large regions of the SUSY parameter space, which we 
outline in detail. We also give, within this setting, some results on the cosmo- 
phenomenologically preferred stau coannihilation region. Finally, we propose a 
parametrization, in terms of the SUSY input parameters, of the common SUSY 
mass appearing in the leading log and mass insertion approximation formula for the 
charged lepton flavor violating decay rates, which fits our full renormalization group 
results with high precision. 
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1 Introduction 



In the Standard Model (SM) with left-handed massive neutrinos, lepton flavor violating 
(LFV) processes, such as the muon decay into electron and photon, \i — > ej, or tau 
decay into muon and photon, r — > /i7, are extremely suppressed: the branching ratios 
for these processes are typically less than 10~ 54 pQ, henceforth far beyond experimental 
reach. However, this situation changes drastically if the SM is embedded in its (minimal) 
supersymmetric extension, including heavy right-handed neutrinos, responsible for the 
generation of low energy neutrino masses via the see-saw mechanism (MSSMRN) [21 E]- 

The compelling evidences for neutrino mixing, coming from solar and atmospheric neu- 
trino, and KamLAND experiments [I], indicate that the lepton charges are not conserved 
quantities: as a consequence, it is not possible to choose a flavor and gauge interaction 
diagonal basis in the lepton sector (see, e.g., [5J). 

The off-diagonal elements of the supersymmetry soft-breaking (SSB) terms in the slep- 
ton mass matrices, (m~)^., (mf)^., and trilinear scalar couplings, (A e )^. {i ^ j) represent 
new direct sources of lepton flavor (charge) violation which are not suppressed by the 
extremely small light neutrino masses. The resulting LFV decay branching ratios, though 
highly dependent on the structure of the theory (e.g., the right-handed neutrino masses, 
neutrino Yukawa couplings, slepton soft SUSY breaking terms etc.) are typically in the 
range of sensitivity of the current and planned experiments searching for such decays. 

A critical issue in the study of the predictions for the LFV decay rates in the SUSY 
theories incorporating the see-saw mechanism of neutrino mass generations is the relation 
between the high energy input parameters and the low energy phenomena. This relation is 
mainly dictated by radiative corrections encompassed in the renormalization group (RG) 
equations. 

A common approach to RG evolution in the context of LFV is to resort to the simple 
leading logarithmic approximation. In [3] it has been shown that, in the specific case of 
sequential dominance in the neutrino Yukawa couplings, the leading log approximation 
may lead to an under evaluation of the off-diagonal entries of the slepton mass matrix. 
The aim of the present work is to show in a general setting what are the effects of 
taking into account the full RG running and to provide evidence of what is the parameter 
space where the leading log approximation ceases to give a reasonable approximation to 
the full result. We do not make any particular assumption about the neutrino Yukawa 
coupling matrix, but instead we resort to a best fit procedure, which allows us to sit 
on a (particular) viable solution, compatible with the known parameters of the neutrino 
sector. We also calculate, within this RG approach, the absolute value of some LFV 
decay rates in a cosmologically and phenomenologically allowed parameter space region 
of a minimal supergravity (mSUGRA) MSSM, namely, the one where neutralino-stau 
coannihilation processes take place. Finally, we present an effective parametrization of 
the common SUSY mass m s appearing in the commonly used approximate formula for 
the lj — > la [rai j > rriij decay branching ration, Z 3 = r, h = h = e, in terms of 
the high energy input parameters of mSUGRA, the universal gaugino mass Mx/% and the 
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common SSB scalar mass mo, both taken at the GUT scale. We demonstrate that our 
parametrization allows to reproduce with high precision the exact RG result in a large 
region of relevant parameter space. 

The paper is organized as follows: in the next Section, we discuss the current phe- 
nomenological constraints on the SUSY parameter space, on neutrino masses and mixings 
and on LFV processes. In Section El we illustrate our numerical approach and methodol- 
ogy, while in Section |H we discuss the relevant RG equations. Our results are presented 
in Section El Finally, Section El is devoted to conclusions. Details regarding the formalism 
used for the calculations performed in the present study are given in Appendices 1X1 and iBl 

2 Phenomenological Setting 

The main ingredients entailing non negligible rates for LFV processes in the context of the 
MSSMRN are (1) the structure of the SSB terms and (2) the neutrino mixing in the lepton 
sector of the theory. We devote this section to the discussion of the SUSY breaking sector 
(Sec. EH}, which is currently restricted by both cosmological and accelerator data 0, 
and the problem of determination of the high energy structure of the Yukawa couplings 
of the theory (Sec. I2.2JI . It is well known that low energy data allow to reconstruct only 
partially the high energy couplings of the theory. Finally, we will briefly discuss how LFV 
processes arise in the context of MSSMRN theories (Sec. 12. 3|) . 

2.1 SUSY Parameter Space 

In the MSSM, the Standard Model fields are promoted to superfields and the theory is 
defined by a superpotential W and by a soft SUSY breaking Lagrangian j£? so ft- This last 
part results after a formal integration over the hidden sector fields which break SUSY. 
In its most general setup, Jz^ ft depends a priori on more than one hundred parameters, 
which are, however, constrained both by theoretical arguments (e.g., naturalness in the 
Higgs sector) and experimental data (e.g., absence of flavor changing neutral currents). 
The specific form of J*f SSB relies on the mechanism which communicates SUSY breaking 
from the hidden to the visible sector. If one resorts to non-trivial assumptions about 
high energy physics, such as Grand Unification or specific supergravity Kahler potentials, 
a further reduction of the number of parameters can occur. One of the most common 
settings is the so-called minimal supergravity (mSUGRA) 7\. The auxiliary chiral field 
giving mass to the gauginos is supposed to lie, in mSUGRA, in the trivial representation 
of the symmetric product of two adjoints of the GUT gauge symmetry group, hence 
generating a universal gaugino mass, M\/i. Moreover, the scalar soft breaking part of 
the Lagrangian depends only on a common scalar mass mo and trilinear coupling A , as 
well as on tan (3. Fixing an extra sign-ambiguity in the Higgs mass term fi completes the 
parameter space of mSUGRA, which then reads 

Mi/2, rn Q , A , tan/5, sign(/i) . (1) 
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These parameters are restricted by accelerator and cosmological constraints. An obvi- 
ous accelerator constraint comes from negative searches for superpartners at LEP2 [S]. In 
mSUGRA models, however, the most stringent bounds come from indirect phenomenolog- 
ical implications of SUSY. For instance, the limits on the almost SM like lightest CP-even 
neutral Higgs boson put a strong constraint at > 114.1 GeV, which, being m^(m x ) an 
increasing function of its arguments, translates into a lower bound on m x jU] . Analogously, 
the inclusive branching ratio BR(6 — ► 37) ID] receives SUSY contributions proportional 
to tan/3 and inversely proportional to m x . The current experimental data [TJ], combined 
with the SM theoretical uncertainties, give the overall bound 

2.16 x 1(T 4 < BR(b -> s 7 ) < 4.34 x 1(T 4 . (2) 

The bound of Eq. (j2J) translates into a lower bound for the mass of the LSP. We also 
mention the constraint coming from the deviation 5a M of the measured muon anomalous 
magnetic moment from its SM value in this case the current theoretical uncertainties 
in the SM computations, mainly due to the evaluation of the hadronic vacuum polarization 
contribution, make it difficult to draw a bound from this quantity. Therefore, this bound 
is not taken as a constraint in the present work. 

Since supersymmetric models with conserved .R-parity offer a natural (i.e. stable, 
massive and weakly interacting) candidate for the non-baryonic matter content of the 
Universe (Dark matter), namely the lightest neutralino \i ft is natural to require that the 
relic density fraction Q x of neutralinos falls into the cosmologically preferred range. In 
this respect, the recent data from WMAP ^3] greatly increased the accuracy of previous 
determination of the Dark Matter content of the Universe, indicating that 

o h 2 - n ii2fi +0 - 0081 cvi 

"CDM 11 — U.-L-LZ,U_ 0.0091 • \° ) 

As the lower limit can be evaded under the hypothesis of the existence of another cold dark 
matter component besides neutralinos, we take here as a constraint only the upper bound 
fl x h 2 < 0.13. Imposing on the mSUGRA parameter space (P) the constraints coming both 
from cosmology and from accelerator experiments typically reduces the viable values of 
the parameters m and M1/2 to very narrow strips [To] . 

The major problem, in the framework of mSUGRA, is that the annihilations of bino- 
like neutralinos is generally not sufficiently efficient, and the resulting Q x exceeds the 
previously stated upper bound. Specific relic density suppression mechanisms are needed 
in order to achieve reasonable values for Q x . Exhaustive investigations of the mSUGRA 
parameter space showed that only four regions of the parameter space are compatible 
with Eq. (jHJ) [12]. (1) A bulk region where the neutralino is sufficiently light and no 
specific suppression mechanism is needed; this region is, however, severely restricted by 
increasingly accurate accelerator bounds. (2) A coannihilation region, where the light- 
est supersymmetric particle (LSP) is quasi degenerate with the next-to-LSP (NLSP). In 
this case, the freeze out of the two species occurs at close cosmic temperatures, and co- 
annihilations between the NLSP and the LSP can drastically reduce the relic density of 
the latter. (3) A focus region, at high m , close to the region excluded by the absence of 
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radiative electroweak symmetry breaking, where the LSP gets a non negligible fraction 
of Higgsino, thus enhancing direct annihilation into gauge bosons. (4) A funnel region, 
where 2m x ~ m^, tua being the mass of the CP-odd neutral Higgs boson of the MSSM. 
In this region, which opens only at large tan f3, resonance effects suppress the relic density 
through direct s-channel pole annihilations. 

In this paper we will only exploit the coannihilation region, since the coannihilation 
mechanism can occur at every tan/3, and also takes place in a relatively "stable" SUSY 
parameter space region, contrary to the "focus" and "funnel" regions. Moreover, it has 
been recently shown ^B] that in supergravity models where some relation exists between 
the trilinear and bilinear soft supersymmetry breaking parameters Aq and Bq {e.g., the 
Polonyi model or particular no-scale models) only the coannihilation region survives after 
the cosmological and phenomenological constraints are applied. 

In general, at a fixed value of tan/3, the stau coannihilation strip lies at values of 
m ~ aMi/2 + b. The minimum Mi/ 2 value reflects the lower bound on m x , dictated by 
accelerator constraints, while the maximal iWi/ 2 value is given by the saturation of the 
bound © at m x ~ m NLSP . 



2.2 Neutrino Sector 

The right-handed neutrino Majorana mass term, Mr, and the neutrino Yukawa couplings, 
Y v , produce an effective Majorana mass matrix (see Sec.EJ) for the left-handed neutrinos. 
This is the well-known see-saw mechanism [17 : 

M eff ~ M J u M~ l M v 

~ vl sin 2 PYjM^Yv , (4) 

where v u is vacuum expectation value (VEV) of up-type Higgs field, v u = 246/ v^2 GeV. In 
what follows we do not consider contributions other than that in Eq. (j3J) to the left-handed 
Majorana mass term. 

The matrix M e g is diagonalized by a single unitary matrix according to 

M es S = Uj s M eS U e{i = diag (mi, m 2 , m 3 ) . (5) 

Let us define by U e the unitary matrix which diagonalizes M\M e , where M e is the charged 
lepton mass matrix, 

UlMlM e U e = diag (mlml,m 2 T ) , (6) 

M e = Y e Vd cos P, Y e and Vd being the charged lepton Yukawa couplings and the down-type 
Higgs VEV. Then, the MNSP neutrino mixing matrix has the form: 

U M NSP = UlU e ff 

(C13C12 C13S12 s 13 e~ lS \ 

-C23S12 - S23Si 3 ci 2 e l<5 C23C12 - s 2 3S 13 s 12 e lS s 23 ci 3 diag(e 10 , e l</3 , 1) , (7) 

S23S12 - C23Si 3 Ci2e 1 ' 5 -S23C12 - C23Si 3 Si 2 e l<5 C23C13 / 
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where we have used the standard parametrization of U MNS p and the standard notations 
Sij = sin0y, Cjj = cos Qij. In Eq. fl2J), 5 is the Dirac and <ft an d y 9 are the two Majorana 
CP violation phases [THj . 

The solar, atmospheric and reactor neutrino experiments |3] have shown that ^H] 

Am 2 l2 = Am| ~ 7.32 x KT 5 eV 2 , Am 2 3 = Am 2 tm ~ 2.6 x 1(T 3 eV 2 , 
tan 2 #12 ~ 0.41 , tan 2 #23 ^ 1-0 , sin#i 3 < 0.2. (8) 

This information does not allow us to distinguish between the three possible types of 
neutrino mass spectrum: hierarchical (m 3 3> m2 3> ^i), inverted hierarchical (m 3 <C 
ni2 ~ mi) or quasi-degenerate (m 3 « vri2 ~ mi). We will assume that the mass spec- 
trum of light neutrinos is hierarchical and therefore m 3 ~ a/ Am 2 tm and m2 ~ a/AttIq. 
Additionally, we assume that all mixing angles lie in the interval < #12, #23, #13 < vr/2. 

We summarize next all assumptions we make in order to fix the Yukawa couplings at 
the GUT scale: 

• All Yukawa coupling constants are taken to be real. 

• The right-handed neutrinos are degenerate, Mr = Mr 1. 

• The neutrino mixing angles are in the interval (0, 7r/2). 

• The CHOOZ angle is given by tan# 13 = 0.1; the other parameters are fixed to the 
values given in Eq. (jSJ) 

• The neutrino spectrum is hierarchical: m 3 = a/ Am 2 tm , m 2 = a/ Ami, mi <C m 2 . 

2.3 Lepton Flavor Violation in the MSSMRN 

The see-saw mechanism requires the existence of heavy right-handed neutrinos as well 
as of the neutrino Yukawa coupling Y v . Lepton flavor violation is then unavoidable because 
in general Y e and Y v cannot be diagonalized simultaneously. In the Standard Model with 
the right-handed neutrinos, LFV processes, though allowed, are suppressed well below 
the sensitivity of any planned experiments. In supersymmetric models, the situation is 
different because there is a new possible source of lepton flavor violation - the soft SUSY 
breaking Lagrangian, Jz^ ft> whose lepton part has the form: 

- ^f aoft = (mDijllZj + im 2 ~) ij e f R ^ Rj + (rr4) i:j v* R -v Rj 

+ (^(A e )ijHd£ Ri Lj + (AJijHjfaLj + h.c^j . (9) 

Non-vanishing off-diagonal elements in the slepton mass matrix are a source of lepton 
flavor violation. If they are present they must be relatively small in order to satisfy 
the stringent experimental bounds on the rates of LFV processes jH]: BR(/i — > ej) < 
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1.2 x KT 11 , BR(/i -»• 3e) < 1 x KT 12 and BR(r -> ^7) < 1.1 x KT 6 . In mSUGRA models 
it is assumed that, at the GUT scale, the slepton mass matrices are diagonal and universal 
in flavor, and that the trilinear couplings are proportional to the Yukawa couplings: 

{ml)ij = {ml)ij = (mDij = (m|) - = (mj)^ = (m|) y = (%m 2 , 

m H d = m H u = m l , (10) 
{A v )ij = A^iYy)^ , (A e )ij = AQ(Y e )ij , [A u )ij = Ao(Y u )ij , (Ad)^ = A^iY^)^ . 

However, soft-breaking terms are affected by renormalization via Yukawa and gauge inter- 
actions, so the LFV in the Yukawa couplings will induce LFV in the slepton mass matrix 
at low energy. The RG for the left-handed slepton mass matrix is given by 

^(m 2 7 ).. = fi-fj- (m 2 T ) .. +^-r\(m 2 L YX + YX^l)- 
^dn K L> ^ ^dfi K mssm 16vr 2 L v lv l} ■> 

+2 (Yjrr%Y v + m 2 H YjY v + A[A V ) .. 

where the first term is the standard MSSM term which is lepton flavor conserving, while 
the second one is the source of LFV. In the leading log approximation and with universal 
boundary conditions, Eq. (jl()j) . the off-diagonal elements of the left-handed slepton mass 
matrix at low energy are given by 

These off-diagonal mass terms generate new contributions in the amplitudes of LFV pro- 
cesses such as /i — > e + 7 and r — > /x + 7. 



3 Method of Numerical Computation 

In order to perform a calculation of the rates of the lepton flavor violation processes, we 
have to use Yukawa coupling constants that reproduce correctly the low energy fermion 
masses and mixings. As we have already discussed in Sec. 12.21 we have assumed that 
all Yukawa coupling constants are real, and that the down-type quark, Y d , the charged 
lepton, Y e , and the right-handed neutrino mass matrices, Mr, are diagonal in the flavor 
basis at the GUT scale. This implies that up-type quark Yukawa matrix is not diagonal 
and is the source of the CKM mixing matrix: 

\ 

x C • (12) 

In the standard parametrization the CKM matrix reads: 

(C12C13 S12C13 
-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 I , (13) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 
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Experimental values 


m u = 


2.80 MeV 


m c = 0.60 GeV 


m t = 166.71 GeV 


m d = 


4.75 MeV 


m s = 90.50 MeV 


m b = 2.90 GeV 


m e = 


0.511 MeV 


m M = 105.0 MeV 


m T = 1.746 GeV 




0.22 


V cb = 0.041 


V ub = 0.0035 


Am 


= 7.32 x 10~ 5 eV^ 


Am 2 tm = 2.6 x 10" 3 eV^ 




tan 2 6< 


= 0.41 


tan 2 atm = 1.0 


tan 2 ^13 = 0.01 



Table 1: Conventional experimental data at Mz- 

where = cos% and Sj-,- = sin 6V,-, i,j = 1,2,3 being generation indices, and we have 
neglected the CP violation phase. 

Let us note that the assumption that the matrix of the charged lepton Yukawa cou- 
plings, Y e , and the right-handed Majorana mass matrix, Mr, can always be simultaneously 
diagonalized at the GUT scale is not fulfilled in all GUT theories. Assuming that Y e and 
Mr are diagonal implies that the solar and atmospheric neutrino mixing angles are gen- 
erated essentially by the neutrino Yukawa couplings, Y u , since Y e and M R do not change 
substantially by the RG effects. In other words, the MNSP matrix is essentially the diag- 
onalizing matrix of the see-saw light neutrino mass matrix. If the charge lepton Yukawa 
matrix is not diagonal, its diagonalization will contribute significantly to the neutrino 
mixing in the weak charged lepton current. 

Under the assumption made of real Y u , the latter depends on 9 parameters at the 
GUT scale. In the numerical computation we calculate the quark, charged lepton, the 
two heavier neutrino masses, and the CKM and MNSP mixing angles, using a set of 21 
free parameters contained in Y u (Eq. Y d , Y e , Y v (6 + 3 + 3 + 9 = 21) at the GUT 

scale. We treat Mr as an input parameter. We use the one-loop renormalization group 
equations 1 to obtain the values of these parameters at the weak scale, set to the mass of 
the Z boson (Mz = 91.188 GeV). In order to find the "best possible fit", we define a 
quantity called b.p.f., which depends on the ratios of the fitted masses and mixing angles 
to the experimentally determined masses and mixing angles at the energy scale of Mz 
(see in Tab. HJ): 

(14) 

Here (/) are the fitted masses and mixing angles with a set of Yukawa coupling constants 
and / exp are the corresponding experimental values. We select Yukawa coupling constants 
which give a minimal value of b.p.f. less than 10 -2 . In other words, the Yukawa coupling 
constants which we will use for the numerical calculation of the lepton flavor violation 
decay rates can fit the 17 low energy fermion masses and mixing angles within an average 

the renormalization group equations for the gauge coupling constants and gaugino mass terms we 
have included a part of the two-loop contributions as well. For further details see Appendix 1X1 



b .p.f. s £ h jfl 

V J cxp / 
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deviation from the experimental values of exp <yi0- 2 /17) « 1.025, i.e., we can reproduce 
the indicated low energy values with a deviation from the measured ones which on average 
is less than 3%. 



4 Renormalization Group Equations from the Uni- 
versal Right-Handed Neutrino Scale 

In this section we will discuss the renormalization group equations treatment of the 
neutrino Yukawa couplings, the right-handed Majorana mass term and of the trilinear 
coupling term, A v . From the gauge coupling unification scale - the GUT scale - to the 
universal Majorana right-handed neutrino mass scale we use the MSSM renormalization 
equations which are given in Appendix El Below the see-saw scale, Mr, the right-handed 
neutrinos are integrated out. Thus, Mr, Y v and A v are not present in our set of renor- 
malization group equations below Mr, since they are no longer physically relevant: their 
effects are in fact encompassed into the running of the effective neutrino mass matrix 
M e ff below that scale. This implies that the RGE for the up-type and leptonic Yukawa 
/3-functions below the universal right-handed Majorana mass scale and up to the super- 
symmetry breaking scale, M susybk , 

M susybk = v /m~m~' , (15) 
where mr and mj are the stop masses, read: 

levr 2 ^ = [-^g!-Sg!-fg! + 3TT(Y u Y^Y Uij 

+3 {Y u Y^Y u ) l3 + {Y u Y ] d Y d ) i3 , (16) 
16tt 2 ^ = |-^-3^ + 3Tr(F d yJ)+Tr(Y e F e t)}y ej . + 3 (Y e Y*Y e ) .. ,{YI) 

16tt 2 |k).. = -(^gimf + egilM^s^-yss.j 

+ {mlYX + YlY e n%).. + 2 [Y^Y, + m 2 H YX + A\A e ).. (18) 
^ 2 J t A ^ = {-^i 2 -3^ + 3Tr(FiF,) + Tr(Y e tY e )|A eij 

+2 j-^Mi + 3g 2 2 M 2 + 3Ti(YjA d ) + Tr (Y^Ae) } Y £i . 

+4(y e y e U e ).. + 5(A e y e ty e ).. , (19) 

l^ 2 j t A Uij = !~g* - 3g 2 2 - ^ g l + 3Tr(YXu)\A Ui . 

+2 ^9{M X + 3g\M 2 + ^gjM 3 + 3Tr(Y^A u ) } Y Uij 
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+4(y tt yjA) ij + b{A u Y^Y u ) tJ + 2{Y u Y}A d ) ij + {A u Y^Y d ) t3 , (20) 
16vr 2 |(mij = -{^gl\M^ + Qgl\M^ + \ g lS 

+6 Tr (m|r u tr u + yftmg + m^)Y u + A{A U ) . (21) 

Note that the other renormalization group equations do not change in this energy 
interval, i.e., we apply the same equations for the gauge coupling constants, gaugino 
mass terms, Y d , m%, m~, m~, m~, A d and m\, as in Appendix 1X1 

In an analogous way, we take into account the running of the effective neutrino mass 
matrix considered as a non-renormalizable term [20] from the universal right-handed Ma- 
jorana neutrino mass scale to the supersymmetry breaking scale, Eq. (JT5J) : 

167r 2 ^M eff = i^—gl - 6gj + 6 Tr (Y^Y u )j M cff + (y e t F e ) T M eff + M eE (Y}Y e ) . (22) 

Below the supersymmetry breaking scale we use the usual renormalization group equa- 
tions for the Standard Model and non-supersymmetric version of the see-saw mass matrix 
renormalization equation (see the set of renormalization group equations in |21j). We 
use as well the relation between the the Weinberg-Salam Higgs field, WS , self-coupling 
constant A and the gauge coupling constants 

\ = \{gl + gl) • (23) 

This relation is a consequence of the underlying supersymmetric structure of the theory 
(see, e.g., [22!). 



5 Results 

We present in this section the main results of this paper. We begin discussing the correc- 
tions to the leading log approximate values of the off-diagonal entries of the slepton mass 
matrix stemming from the full RGE running. We then review the corresponding effects 
on LFV processes, and discuss the role of the sign of the trilinear coupling A and of \i. 
Finally, we give predictions for LFV decay rates, within the present best-fit approach, in 
the coannihilation region of mSUGRA models and present an effective parametrization of 
the common SUSY scale m s appearing in the "short-hand" formula for BR(/j — > lj) (see 
Eq. (|25jl ). We include in our results predictions at low, intermediate and large tan/3. In 
particular, we took tan/3 = 3 in order to illustrate the effects we find in the extremely 
low tan/3 regime, though in some regions of the parameter space this value is ruled out 
by the LEP2 bound on rrih- 
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5.1 The Left-Handed Sleptons Mass Matrix 



We investigate the effect of full RG running in the off-diagonal entries of the slepton mass 
matrix m~ by numerically evaluating the ratio of the exact result and of the leading log 
approximation result which reads 

( m ^ LL % s -T^< 6m « + M «> ( y ' y -)* '°s psf ) (24) 

We study in Fig. [T]the element (2,1) of the quantity \m\{RG) / m\(LL)\ 2 as a function of 
the common right-handed neutrino mass scale Mr for two representative values of tan j3 = 
3.0, 30.0. We first notice that, in all cases, full RG running yields an increase in the off- 
diagonal slepton mass matrix elements with respect to the leading log approximation. 
Secondly, increasing Mr gives rise to larger effects, whose size increase with the high 
energy common gaugino mass M\/2- This is expected, since the effect of Mi/ 2 in the 
running is completely disregarded in the leading log approximation. For Mr up to 10 9 GeV 
we get a maximal increase of one order of magnitude at M\/2 — 1 TeV, while for Mr « 
10 14 GeV the error one makes taking the leading log approximation amounts to two 
orders of magnitude. Finally, a remarkably weak dependence on tan (3 is found, as is also 
indicated by the comparison of the two panels in Fig. ^ 

In Fig. El we study the (3,2) and (2, 1) entries of the same ratio at tan/3 = 10.0 at 
various values of m , setting Mr = 10 11 GeV. We notice that increasing m reduces the 
effects of the full RG running: for the maximal value of mo = 400 GeV we consider, the 
ratio of interest [m~(i?G)/m~(LL)] changes at most by a factor of ~ 2. 

We then conclude that the leading log approximation is not accurate for the evaluation 
of the off-diagonal entries of the slepton mass matrix in the low itlq and large Min 
regime. The degree of accuracy is further worsened with the increasing of the right- 
handed neutrino mass scale Mr. 

5.2 Lepton Flavor Violating Processes 

In this section we turn to the main topic of interest in the present work: the size of the 
corrections induced by the full RG running in the rates of LFV processes. We compare here 
the results obtained inserting the leading log approximation for (m~(LL))^, Eq. (J2U), in 
the full mass matrix formulae P2J given in Appendix El see Eqs. ()70|) - (fTH|) . with the exact 
RG running result for the LFV decay branching ratios. We would like to emphasize that 
the results for the off-diagonal slepton mass terms derived in leading log approximation 
and by exact RG running are included in the slepton mass matrix which after that is 
diagonalized and the mass eigenstate scalar leptons and the corresponding mixing matrix 
determined. The letter are used in the calculations of BR(/j — > ^7). Note that the 
branching ratios of the LFV processes depend on the eigenvalues of the the slepton mass 
matrix and on the elements of the matrix which diagonalize it (see Appendix |B}; they do 
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Figure 1: The ratio of the (2,1) element of m~^(RG) /m'-(LL) 
tan/3 = 3.0 (a) and tan/3 
H > 0. 



as a function of M\/2 at 
30.0 (b). In the two panels we fixed m,Q = 100 GeV, Aq = and 



not depend directly on the off-diagonal entries of the slepton mass matrix. Therefore we 
will in general get results which differ from the ones expected on the basis of the results 
reported in the preceding section. 

In Fig. H2 we plot as functions of Mi/ 2 the full RG results (solid red lines) and the 
results obtained within the leading log approximation (dashed green lines) for two values 
of m , m = 100; 300 GeV, at tan/3 = 30.0 and M R = 10 11 GeV. We see that, in 
agreement with was obtained in the preceding section, the increasing of m reduces the 
RG corrections: the RG flow is mainly driven by m and the corrections due to Mi/ 2 
become less important. We get, for m = 100 GeV, a correction of roughly one order of 
magnitude at large M\/2 ~ 1 TeV, while for the m = 300 GeV case the effect reduces to 
a factor of two. We notice, comparing panels (a) and (b), that the difference between the 
exact RG results and those obtained in the leading log approximation is larger at lower 
values of tan (3. 

Fig. 0] shows the effect on BR(/i — > e^y) of changing the right-handed heavy neutrino 
mass M R as a function of M 1 / 2 - We used m = 100 GeV, A = 0, /i > and tan /3 = 30.0, 
taking M R = 10 9 GeV, 10 12 GeV and 10 14 GeV. We find that the size of the difference 
between the exact RG result and the leading log approximation one depends weakly on 
M R . At large My 2 we always get a discrepancy between leading log and full RG results 
of about one order of magnitude (notice the different scale with respect to the preceding 
figure). Here, again, the effect is enhanced at lower tan/3 (see panel (6)). 
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Figure 2: The ratio of the (3,2) (a) element and of the of (2,1) (b) element of 
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m~{RG)/m~{LL) 
M R = 10 11 GeV, A 



as a function of M\ji at tan/3 
= and u > 0. 



10.0. In the two panels we fixed 



Finally, we would like to emphasize that although the change in the LFV decay rates 
we find using the exact RG results in comparison with those obtained in the leading 
log approximation are within one order of magnitude, it would be mandatory to use the 
former in order to draw conclusions concerning the SUSY sector if any LF violating decay 
will be observed and its branching ratio measured. As it follows from Figs. HJ if we fix 
BR(yU — > ej), for instance, we could make an error of up to ~ 50% using the leading log 
approximation in deriving a lower bound on Mi/ 2 . 



5.3 Exact RG Evolution Effects from \i and Aq 

In this section we deal with two other high energy input "parameters" which are neglected 
in the leading log approximation, namely the sign of the scalar trilinear coupling Aq and 
of fi. We depict in the left panel of Fig. UJ the results we get, in the full RG running 
computation, for A = (solid red line), A = ±1 TeV (resp. lower and upper dashed 
green lines) and A = ±2 TeV (resp. lower and upper dotted blue lines). In the figure 
we fixed for definiteness tan/3 = 10.0, M R = 10 11 GeV and m = 200 GeV. We see that 
instead negative values of the trilinear coupling generically enhance, by up to roughly 
a factor of 5, the LFV branching ratios. Flipping the sign of fi gives rise to smaller 
corrections (panel {b)). In particular, we find that at low m , values of fi > slightly 
suppress the LFV branching ratios, while at larger m a mild enhancement takes place. 
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(a) (b) 

Figure 3: The branching ratio for the LFV process BR(/i — > ej) as a function of M1/2 at 
m.Q = 100 GeV and mo = 100 GeV in the full RG (solid red lines) and leading logarithmic 
approximation (dashed green lines). We set in the figures Mr = 10 11 GeV, Aq = 0, fi > and 
tan/3 = 30.0 (a) and tan/3 = 3.0 (b). 




(a) (b) 

Figure 4: The branching ratio for the LFV process BR(/U — > ej) as a function of M 1 / 2 oi 
Mr = 10 9 GeV, 10 12 GeV and 10 14 GeV in the full RG (solid red lines) and leading logarithmic 
approximation (dashed green lines). We set in the figures tuq = 100 GeV, Aq = 0, /x > and 
tan/3 = 30.0 (a) and tan/3 = 3.0 (b). 
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Figure 5: (a): BR(/i ei) as a function of M 1/2 at A = 0, ±1 TeV, ±2 TeV for tan (3= 10.0, 
M R = 10 11 GeV and m = 200 GeV. (b): The same at m = 100 GeV (lower lines) and 
rriQ = 500 GeV (upper lines) at \i > (dashed green lines) and fj, < (solid red lines), for 
tan/? = 3.0, M R = 10 12 GeV. 



5.4 The Coannihilation region 

We provide in this section some results concerning the branching ratios of LFV in a 
particular region of the mSUGRA parameter space, namely the one where coannihilations 
between the lightest neutralino and the lightest stau (the next-to-lightest SUSY particle) 
concur in reducing the neutralino relic density within the current observational cold dark 
matter content of the universe. In Fig. |)]we study the BR(/i — > ej) (a) and BR(r — > yrf) 
(b) as a function of the neutralino mass m x along the "coannihilation corridor" at fixed 
tan/3 = 30.0 for three different values of M R = 10 9 GeV, 10 11 GeV and 10 13 GeV. 
We recall that in mSUGRA m x §AMi/ 2 - Notice that for sufficiently low values of 
Mr the computed BR(yU — > ej) is always below the current experimental bounds, while 
putative lower bounds on m x can be drawn for larger Mr. In the case of BR(r — > /ry) we 
always get results far below (at least two orders of magnitude) the current experimental 
bound. As was natural to expect, and as is shown in Fig. [7| (a), lowering tan/? yields a 
quadratic suppression in the LFV Branching ratios, which appear to be in all cases well 
below the current experimental sensitivity. Finally, in Fig. O (&) we summarize our results 
showing the Mu2 ranges dictated by cosmological and phenomenological requirements. 
These constraints, leading to smaller allowed M\/2 for fi > indicate that this last case 
is favored to produce larger LFV decay rates, which may be in the range of sensitivity of 
the next generation experiments. 
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Figure 6: ( a): BR(// — ► ey) as a function of m x in the stau coannihilation strip at tan (3 = 30.0, 
/i > andA = 0, for M R = 10 9 GeV, 10 11 GeV and 10 13 GeV. (b): the same for BR(r -» pry). 




150 200 250 300 350 400 300 400 500 600 700 800 900 1000 1100 

m x < GeV > M 1/2 (GeV) 

(a) (b) 

Figure 7: ( a): Same as in Fig\$ (a) but at tan/3 = 3.0. (b) BR(/u — ► cy) along the cosmologi- 
cally and phenomenologically viable coannihilation strips for both signs of p and various tan/3. 
Results at tan (3 = 3.0 are ruled out by the bound on set by LEP2, but we nonetheless report 
them for completeness. 
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5.5 A Candidate for the Effective SUSY Mass m s 



The branching ratio for the processes lj — > la is often quoted in the literature in the mass 
insertion and leading log approximations as 



a 3 



BR(/j - la) ^ g-^— tan^ (3 , (25) 

where m s is a typical mass of superparticles, a ~ 1/128 and (m~(LL)) .. is given by 
Eq. (|24|). The problem with this formula is that there is no prescription for m s in terms 
of the fundamental SUSY and soft breaking parameters of the theory, and the dependence 
on m s is so strong that the predictions for BR(/j — > la) depend drastically on what one 
uses for m s : m s = m , or m s = m„, or m s = JlMy^. The three indicated choices will 
give completely different predictions for the branching ratios of interest. Notice that if 
we underestimate m s by just a factor of two we are overestimating the branching ratios 
by more than two orders of magnitude. 

We have found that the expression (|27)jl with (m~(LL)) given by Eq. (J24|) represents 
an excellent approximation to the exact RG result if for m s one uses 

ml ~ 0.5 ml My 2 {ml + 0.6 M 1/2 ) 2 . (26) 

In Fig. IS] we compare the predictions for m| obtained from (|25|) and the exact RG 
result for BR(/j — > ^7) with that given by Eq. (|26|). In general, we find that our fit formula 
(J2HJ) slightly underestimates m s at low M1/2 and large mo, and somewhat overestimates it 
at large My^ and small mo. The deviations are, however, always relatively small, and the 
overall dependence of m| on mo and Mi/ 2 is everywhere reproduced rather accurately. 



6 Conclusions 

We have shown that the off-diagonal elements of the slepton mass matrix obtained by 
RG running deviate significantly - at low mo and large M1/2 - from the values given 
by the leading logarithmic approximation. The exact result is always larger than the 
approximated one and this difference increases for larger Mr. As a consequence, in the 
leading log approximation the predictions for \i — > and r — > ^7 can be underestimated 
by a factor ~ 10. We pointed out that, even when this factor is smaller, it is relevant if 
one wants to use BR(/i — > e-y) and BR(r — ► ^7) to constrain the SUSY parameter space. 
These branching ratios were found to depend on the sign of Aq and of /1 as well. 

We also studied the behavior of BR(/x — > ej) and BR(r — > /17) along the cosmologically 
determined "coannihilation strips" . Our results show that the present experimental bound 
on /i — > e7 can exclude only regions with large Mr (Mr > 10 13 GeV), while r — >• ^7 is 
always allowed. Finally, we proposed a simple parametrization in terms of m and My 2 
for the effective SUSY mass parameter m s which enters into the simple expression for 
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M 1/2 (GeV) 

Figure 8: A fit to m§. The lines correspond to Eq. \2b]) . while the points to the value extracted 
from the full RGE numerical results and Eq. \25\) . In this plot we use Mr = 10 14 GeV, tan/3 = 3 
and Aq = 0. 



BR(7j — > lij) derived in the leading log and mass insertion approximations (see Eqs. (124)1 - 
(J2EJ))- This allows to reproduce the exact RG results for BR(/j — > ^7) with high precision. 
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A Renormalization Group Equations 

From the GUT scale M GUT down to the universal M R scale we use the one-loop MSSM 
renormalization group equation [2I| 2 of the gauge coupling constants, the three gaugino 
mass parameters, the Yukawa coupling constant matrices, Y u , Yd, Y e , Y u , and the right- 
handed neutrino mass matrix, M R , and the soft supersymmetry breaking mass terms and 
trilinear parameters, A. 

The gauge coupling constants (especially gi running effect) and the gaugino mass 
terms (as well M 2 running effect), play important role in the calculation of the lepton 

2 We adopt S. P. Martin and M. T. Vaughn convention , «-e., the sign of the terms proportional to 
the gaugino masses in the trilinear parameters, A, are different from those used in |23j . 
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flavor violation processes, therefore a part of the two-loop (3 functions are adopted in our 
computations. We present at first the non-supersymetric part of the renormalization from 
the GUT scale M GUT to the universal M R scale: 

1R 2d 33 3 gf /199 2 27 2 88 2 \ 

1671 37 9i = ~r 9i + \ -*T9i + ~r92 + ~r93 ] > ( 27 ) 



di " 5 " x ' 16tt 2 V 2V' 5 ^ 5 

^ ^ - g 2 + -^2 

Jt 93 - ~ 393 + 16^\5 



i^±g 2 = gl + S-JU + Kgl + Mgi) , (28) 



16tt4 ^3 = -3 5f + ( \g\ + 9^ + Ugj ) , (29) 



16vr 2 |F U!J = [-^-Sgl-fgl + STriY^ + TriY^Y^ 

+3 (YXYuh + (YuYjY^ , (30) 

le- 2 ^. = {-^-3, 2 2 -^ 3 2 + 3T^ 

+ 3 (Y d YjY d )ij + (Y d Y*Y u )ij , (31) 

i6vr 2 |F ei3 = {-^ 2 -3^ 2 + 3Tr(y d yi)+Tr(y e y e t)}y eij 

+3 (y e y e ty e ).. + (y e y i ty,).. , (32) 

16^ 2 |n !3 = {-^ 1 2 -3^ 2 + 3Tr(r u r u t) +Tr ( nFj t)|y ;i3 

+3 (ni^n) y + (y^y e ) y , (33) 

^ 2 jM RlJ = 2(M R YX) ij + 2(Y»YjM R ). j . (34) 

Here are the renormalization group equations of the gaugino mass terms, the soft mass 
terms and the trilinear parameters as well the up-type and down-type Higgs mass terms, 
respectively. Note that we take into account the two-loop (3 function accuracy only for 
the gaugino mass terms: 

,d w 66 2g\ f 199 2 , 

167T 2 _Mi = — ( 7l 2 M 1 + -^-^— 5 2 (2M0 



5 ai 16tt 2 5 



+^ 2 (M x + M 2 ) + ^ 3 2 (Mi + M 3 ) \ , (35) 



-2 ^ s , _ o 2,^ , 2 02 f 9 2 



16tt 2 -M 2 = 2g 2 2 M 2 + -^^-gi(M 1 + M 2 ) 

+253 2 (2M 2 ) + 24( ?3 2 (M 2 + M 3 )} , (36) 
16^4 M 3 = -6^M3 + -H(^ 2 (M 1 + M3) 



° ^ ° 16tt 2 [ 5 

9( ?2 2 (M 2 + M 3 ) + 14( ? 3 2 (2M3)} , (37) 
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16n 2 ^A e .. 
dt lJ 



dt 13 



16 ^K).- = -(^i 2 |M 1 | 2 + 6, 2 2 |M 2 | 2 + |, 3 2 |M 3 | 2 )^ + ^ 1 2 5^ 

\ j ij 

+2 {Y^mlY u + m 2 H Y^Y u + 4A U ) . . 

+2 (yjm|y d + m^yjy, + a\a a ) , 

+2 (mlY u Y^ + Y u Y^ml).. 
+A(Y u mlY^ + m 2 H Y u Y^ + A u Al) ij , 

= -(^ 1 2 |M 1 | 2 + f^|M 3 | 2 )^ + ^ 1 2 5^ 

+2 (m|y,yi + y d yX~).. 

+4 (y d m§yj + m^y.y] + A d A\) , 

+ [m\YlY e + mlyjn + y^m? + Y}Y v rr%).. 
+2 (r e tm|r e + m 2 H YX + A\A e ).. 
+2(Y}m?Y u + m 2 H YX + AiA l/ ) ij , 

+4(F e m 2 F e t + m 2 , d y;F e t + A e At) y . , 
16tt 2 ^ (m 2 ,) .. = 2(r4Y v Y^ + Y v Y^r4).. + A(Y v mlY^ + m%^ 



—fx - zgl + 3Tr(r;r d ) + Tr(y e ty e ) } a 6ij 

+2 j^Mi + 3g 2 2 M 2 + 3Tr(Y}A d ) + Tr(y e U e ) \ Y eij 



+4 (n^) y + 5 (A e YX) t] + 2 (y e ^) y + (A e YX) .. , 
+2 | ^ 2 M X + 3£ 2 M 2 + 3Tr(yJ A u ) + TrtfA,) } ^ . 

+4(y„y I M„) ij + 5(A u YXu)ij + 2(y I/ y e U e ) iJ - + (A^y)*, , 
= { -§5i 2 - 3s 2 - y ^ 2 + 3Tr(yJy u ) + Tr(yjy ) } a,, 
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+2 |^<7?M a + 3g 2 2 M 2 + ^ 3 2 M 3 + 3Tr(Y^A u ) + Tr(Y]A v ) } Y Uij 
+A(Y u Y^A u ) ij + h{A u Y,lY u ) t3 + 2{Y u YjA d ) ij + {A u Y ] d Y d ) tJ , (46) 
16vr 2 |^. = j-l^ 2 -3s 2 2 -^ 2 + 3Tr^ 

+2 {^i 2 M! + 3^ 2 2 M 2 + ^ 3 2 M 3 + 3Tr(yjA d ) + TV(y e U e ) | Y dj . 

+4(y d i^Ai)« + 5(^y d V d )^ + 2(Y A YiK)a + (^Xk , (47) 
16 ^|«J = -fe 2 |M 1 | 2 + 6 52 2 |M 2 | 2 )+^ 



+6 Tr (m|Yjy u + l^(m| + m^jy + 4^ 
+2TV(m|y i ty !/ + yt(m| + m 2 ? jy i , + 4^) , (48) 

+6 Tr (m 2 ~YjY d + yj(m§ + m\ d )Y d + 

+2 Tr {m\jX + l?(ro§ + rn^jy + A\A e ) , (49) 



where 



S = Tr(m~ + ra|- 2m| - m| + m§) - m 2 ^ + m 2 ^ , (50) 



and we have used the GUT convention for the U(l) gauge coupling constant, gx, and 
t = ln/i where fi is denoted as the renormalization point. 



B Notations in the MSSM 

In this appendix we set the notations and conventions for the computation of LFV 
processes. Note that we follow the Appendix B of ref. |2Hj . 

The chargino mass matrix has the following form: 



-S? m =[Wn # 2 ~)M C ( ^ )+h.c, (51) 



where 

Mc= ( M 2 V2m w cosP\ 
\ yzrriw sin p \i ) 

This matrix can be diagonalized by two 2x2 real orthogonal matrices L and R according 
to: 

R M c O J L = diagonal . (53) 
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If we define 



XlL 

X2l 



T 



H 



lL 



Xl_R 

X2r 



= o R 



H, 



2r 



then 



Xa X/ 



X 



AR 



(54) 



(55) 



forms a Dirac fermion with mass M~- . 

X-A 



The mass matrix of the neutralino sector is given by 

/ B L \ 



- J£ 



B L W° L H° 1L H° 2L ) M, 



IN 



wl 

\ H 2L ) 



+ h.c. , 



(56) 



where 



M K 



( 



Mi 






M 2 



—mz sin 9w cos (3 mz sin 6w sin (3 \ 
mz cos 6w cos (3 —mz cos 6w sin (3 



—mz sin 9w cos (3 mz cos 8w cos (3 
\ m z sin 6 W sin (3 —m z cos 6w sin (3 





-fi 



-n 







We can diagonalize M N with a real orthogonal matrix h 

N M N T N = diagonal . 
The mass eigenstates are given by 



~o 

Xbl 



(O n ) bc X° cl (S,C = 1,...,4) 



We have four Majorana spinors 



Xb Xbl Xbr i 



x cl = (B L , W L , H lL , H 2L 



(B = l,- -.,4) 



/ 

(57) 

(58) 
(59) 
(60) 



with masses M^o . 

X-B 



The slepton mass matrix can be cast in the following form: 

— — (el , 



2 2+ 



» 2 « 



with 



( m L) 4J 



( m l)ij + m e t ^i ~ m l cos 2/9 ^ + sin 2 0^ , 

(%).. + m\_bij — m| cos 2/5 sin 2 0v^ij , 
f cos /3 



v/2 



(61) 

(62) 
(63) 
(64) 
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We diagonalize the slepton mass matrix, M\ , by a 6 x 6 real orthogonal matrix U l as 



(65) 



U l M 2 L U lT = diag. (m| 

A mass eigenstate is then written as 

Ty = U l Y J Li + Ul ii+3 T Ri , (Y = 1, • • • , 6) . (66) 
Since there are no right-handed sneutrinos at low energies, the sneutrino mass matrix is 



{Mi)ij = (ml) + -m% cos 2/3% . 



ij 2 

The diagonalization is carried out by a 3 x 3 orthogonal matrix U u 



and the mass eigenstates read 



U v MlU vJ = diag. 



ux = U u Xi D Lil X = l,2,3. 



(67) 

(68) 
(69) 



The amplitudes A L R for the processes lj — > 1^ are given by the sum of the chargino 
and neutralino contributions: 



4(c) , An) 



(70) 



—y— 

9-n- 2 Z-^/ m?. 



32vr 2 ^ ml 

A,X v x l 



r L(l) r L(l)* 1 

3AX lAX 6(1 - x AX y 



. n L(l) n R(l)* Xa 



32vr 2 ^ m? 

BY l Y 



TTlj (1 - X AX ) 3 



[-3 + 4x AX - x 2 AX - 2 In x AX ) 



(71) 



1 



-4(1 - 6y By + 3y| y + 2y^ Y - 6y| y \ny BY ) 



, N L(l) N R{l)*^L 

WW m, (l-i/ sy ) 3 ' 



,(c,n) _ ^(c,n)i 



where the dimensionless parameters x AX and y By are defined as 



x 



ML 

Xa_ 

AX — o ' 

m~ 



Vby 



Xb_ 

2 ' 

w 



and the coefficients in Eqs. (|7T|) and (J72|) are given by 



</2 



c, 



H(J) 



-92(0 R ) Al U x4 



(0 L ) A2 U X , 



N 



L(l) 



N 



H(I) 



#2 



<j [-(0„) fla - (O iV ) sl tan^]^ il 



mvy cos (5 



{o N ) m u. 



1 

Y,i+3 



(72) 
(73) 

(74) 

(75) 
(76) 
(77) 

• (78) 
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Using the above mentioned amplitudes, the branching ratio of the decay lj — > ^7 is given 
by 

BR{l -,l ll ) = ^^{\A L \ 2 + \A E \ 2 ) , (79) 

F 

where G P is the Fermi constant and a = e 2 /Ait. 
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